Space-Time geometry and thermodynamic properties of a self-gravitating ball of fluid 

in phase transition 
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A numerical solution of Einstein field equations for a spherical symmetric and stationary system of 
identical and auto-gravitating particles in phase transition is presented. The fiuid possess a perfect 
fluid energy momentum tensor, and the internal interactions of the system are represented by a van 
der Walls like equation of state able to describe a first order phase transition of the type gas-liquid. 
We find that the space-time curvature, the radial component of the metric, and the pressure and 
density show discontinuities in their radial derivatives in the phase coexistence region. This region is 
found to be a spherical surface concentric with the star and the system can be thought as a foliation 
of acronal, concentric and isobaric surfaces in which the coexistence of phases occurs in only one of 
these surfaces. This kind of system can be used to represent a star with a high energy density core 
and low energy density mantle in hydrodynamic equilibrium. 



I. INTRODUCTION 



In relativistic astrophysics and pure general relativ- 
ity the concept of "phase transition" is usually associ- 
ated with very different scenarios. In general relativity 
we have the spherically symmetric collapse of a massless 
scalar field or Yang-Mills field that leads to a formation 
of a black hole the astrophysical context we 

have first order phase transition in neutron star models, 
where we found matter condensation from neutron mat- 
ter to pion-condensed matter and quark matter 0, 0, @ . 
Also we can have a phase transitions in the core of a ro- 
tating neutron star |g. In the last few years, the phase 
transition phenomenon have become important in rela- 
tivistic astrophysics since we can have emission of grav- 
itational waves 0, 111 associated to this transition. Ob- 
servational evidence of this phenomenon can be attained 
with the present Laser Interferometer Gravitational- 
wave Observatory (LIGO). The LIGO-I maximum dis- 
tance for detection of gravitational waves is about 6.4 
Mpc, well beyond the Andromeda galaxy (M31), whereas 
the next generation LIGO-II detectors will probably see 
phase transition events at distances two times longer Q . 

The model studied in this work is indirectly related 
to the topics mention above. Our aim is to study the 
behavior of the thermodynamic properties and the space- 
time structure in a simple albeit important astrophysical 
object, a self-gravitating ball of fluid performing a generic 
first order phase transition. In particular, the Riemann- 
Christoffel curvature tensor is study in some detail for the 
space-time associated the different regions of the matter 
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that presents a first order phase transition. 

To model this system in the context of General Relativ- 
ity, we used the TOV equations (Tolman-Oppenheimer- 
Volkoff) [HllJ- The TOV equations represent the Ein- 
stein field equations for a stationary, spherical symmetric 
system with a perfect fluid energy momentum tensor and 
a local equation of state. 

The most simple and general equation of state to model 
a phase transition in a classical system of particles, is 
the van der Walls equation of state together with the 
Maxwell energy balance construction. The Maxwell en- 
ergy balance construction have been studied in a rigorous 
treatment in systems with two-body interaction potential 
proving the equal-area rule at the transition region in the 
equation of state |0| . Also the non analyticity of the van 
der Walls equation of state together with the Maxwell en- 
ergy balance construction appears to be generic for flrst 
order phase transitions originated in microscopic short 
range forces 0, 0| . 

To solve the TOV equations one needs a smooth equa- 
tion of state. In order to overcame the fact that in a 
first order phase transition we will have a non analytic 
equation of state we shall perform several smoothing in- 
terpolations. Afterward we show that the real situation 
can be obtained as a limit that is independent of the used 
interpolations. 

The article is organized as follows, in section ^1 we 
present the basic equations that model our system, in 
particular we study the TOV equations and present a 
van der Walls like equation of state that models a phase 
transition. In section ITTTl we analyze the special case of 
constant pressure. This case needs to be treated in a 
different form because, strictly, the TOV equations are 
not valid in this case. In section Hvl we present the be- 
havior of the thermodynamic variables, the metric com- 
ponents, the Kretschmann curvature invariant and the 
components of the Riemann-Christoffel curvature tensor 
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inside the star undergoing a phase transition. Finally, in 
section [3 we summarized our results. 



II. STELLAR MODEL FOR A SYSTEM IN 
PHASE TRANSITION 

We consider a static spherical symmetric space-time 
described by the line element 



r^{de^ +si-a^0dip^), (1) 



where v = v{r) and A = A(r). The matter is represented 
by the perfect fluid energy momentum tensor 



P 



Tap = pUaUp + -^{UaUf3 - Qa^) 



(2) 



P = 



{p+ §)im+ ^r^p)c 
r{r — 2m) 



e"^ = 1 —. 



2p' 



(8) 
(9) 
(10) 



The first three equations allow us to solve for the geomet- 
ric mass m(r), the pressure p{r)^ and the energy density 
/o(r), i.e. the thermodynamic quantities of the system. 
The last two equations define the space-time geometry 
through the functions A(r) and u[r) present in the met- 
ric. 

The equation of state ((HJ depends on the particle inter- 
actions of the system. An equation of state that allows a 
phase transition is a van der Walls like equation of state 
of the the form, 



where p is the energy density, is the four-velocity, c 
is the speed of light, and p the pressure measured in the 
local rest frame. We also assume that the pressure and 
the energy density are related by an equation of state of 
the form p ~ p{p). 

The Einstein field equations. 



Rap = -C ( Tap - -j^gapT 



(3) 



where C, — and G is the gravitational constant, can 
be written with the help of and (0) as a set of three 
equations. 
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(11) 



where a and b are characteristic parameters of a partic- 
ular physical system. To solve the TOV equations using 
the equation of state (fTT|l . we have to specify the type 
of particle considered through the parameters a and 6. 
Using the law of the corresponding states is possible to 
obtain an equation of state which is independent of these 
parameters, in other words, valid for any "real gas" . This 
equation of state in the so called reduce quantities (p^, 
Tr, pr) is 
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in which ' denotes differentiation with respect to the co- 
ordinate r. To write the equations above in a simpler 
form we use the function m(r) defined by, 



m(r) = -(1 



(5) 



m(r) = GM{r)/c^ is the geometric mass which has 
dimensions of distance and M{r) is the quantity of mass 
inside a sphere of radius r. With this change of variable, 
we obtain the system of equations known as the TOV 
equations. 



P^p{p), 
m'{r) 



AttG 



-r p. 



(6) 
(7) 



STj-Pr 



(12) 



with Pr, Tr and pr defined by the adimensional relations 



(4) 
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P_ 
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Tr 



T 

tI: 



(13) 



where {pc, Pc, T^) = -Jj, affb) are the values of the 

pressure, energy density and temperature in which the 
first and second derivative of H12|l become zero simulta- 
neously. Also < Pt- < 3. In Fig. ^we show a typical 
van der Walls like isothermal for a system in phase tran- 
sition (gas-liquid) in reduce coordinates. We know that 
the van der Walls like equation of state describes well the 
regions of gas and liquid but in the region of coexistence 
it presents physical inconsistencies, e.g. in the transition 
region we do not have dynamical equilibrium between 
the phases. Therefore, we modify the equation of state 
(|12|1 using the Maxwell energy balance construction |0 . 
The Maxwell energy balance construction allows us to 
calculate the pressure of coexistence po and the coexis- 
tence energy density interval [/Oi,P3] through the system 
of equations. 
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P3 

J Pr (pr) dpr = Po (p3 - Pi) , 
Pi 

(14) 

3pi = PO = ^ 3p3 ■ 

First we fix a temperature value and then we solve the 
system of equations p4f) ( e.g. with the Newton- Raphson 
method for nonlinear systems). Once po is found, wc 
can say that the state equations for our system in phase 
transition between gas and liquid is of the form, 



where f ^ 1 is directly associated with the angle a 
[tana = 2£^/{p^ — Pi — 2x)] of the straight line in the 
transition region. The polynomial coefficients a„ and c„ 
are calculated using the six boundary conditions given by 
the continuity of the pressure, its first derivative, and its 
second derivatives at the junctions. We found that our 
class curve defined in this way is a good approxima- 
tion for the case of a gas-liquid transition. The analysis 
of the case with p = constant will be discussed in section 

EH 

For generality sake, we want that the system of equa- 
tions Ht)llU|l be independent on the characteristic param- 
eters of the gas. We introduce new variables r and F 
define as 



3pr^ < Pr < Pi 



(gas) 



Pr = <^ Po Pi < Pr < P3 (coexistence) . 

- 'ipr'^ Pz< Pr <'i (liquid) 

(15) 

We note that the equation of state H15|) is not of class 
because its derivatives have discontinuities \tl p = p\ 
and p — pz- So, the TOV equations can not be solved 
around these energy densities. Furthermore, in the coex- 
istence region, the equation of state satisfies the equation 
= that brings no physical solutions, e.g. negative 
masses, see equation ©. For this reason, the case of 
constant pressure has to be study independently. The 
problem of smoothness above mentioned can be solved 
using an equation of state of class or smoother in- 
stead of (|15|l . This is done redefining the equation of 
state H15II in the region of coexistence using a virial like 
equation to make, at least, the joints in pi and /03 of class 
. The line Pr = constant in the coexistence region can 
be approximated by a straight line with a small slope an- 
gle (a). In this way, our equation of state for the system 
becomes 



r = Ar, m = AFr 



(18) 



where the quantity A is a scale factor. Furthermore we 
fix our units by choosing 



apc ^ ^ 



1. 



(19) 



Substituting these expressions in equations jnj, and 
© , the internal structure of the star undergoing a phase 
transition is found solving the system of equations 



dpr_ 
dr 



(1-2 rF) 



F 



rpr 
3 



[3/9,. -I- Pr] 



dpr 
dpr 



dF F 
— = - 2—, 
or r 



(20) 



where for each region of the star we have to use the 
isothermal corresponding to the temperature that allows 
a phase transition. 



3pr^ < Pr < Pl-X 

Pi - X < Pr < Pi + X 



5 

E OnP" 
1 

E bnPr Pi + X < Pr < P3 - X 
5 

E CnPr P3-X< Pr < P3+X 
n=0 

L 3-pr 



(16) 



P3 + X < Pr < 3 



where a„, 6„, c„ are constants to be determined, and 
2x is the interval for the coupling polynomials. In our 
case we define the coefficients 6„ so that the straight line 
pass through the point po at the middle of the interval 
[Pi +X,P3-x]: in this way 



bo =Po 



P3 + Pi 
P3 - pi - 2x 



2^ 



P3 - pi - 2x ' 

(17) 



III. ANALYSIS OF THE CASE OF CONSTANT 
PRESSURE 

To perform the integration of the TOV equations 
through the discontinuities of the coexistence region, we 
use the class equation of state (16) with the con- 
stants of the coupling polynomials calculated using the 
six boundary conditions mentioned in section^] The co- 
existence region is modeled, in first approximation, by a 
straight line with a small slope angle (a « 10""^). Then 
we let the angle a go to zero to assure that the abrupt 
change presented in the energy density figure (and other 
figures that depend on the energy density) are indeed ver- 
tical lines. The same result can be obtained in a different 
way. Knowing the pressure and energy densities inter- 
val for the coexistence region (e. g., from the Maxwell 
construction) one can solved the TOV equation indepen- 
dently for the liquid phase as well as for and gas phase. 
Then the end points are joit with a straight line. With 
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this second method we never integrate the TOV equa- 
tions through the discontinuities. Both methods give the 
same results, we use the first one because in our opinion 
the results appear in a more transparent way. 

To determine the behavior of the energy density in the 
coexistence region and the radial interval in which the 
transition occurs in the limit a — > 0, we define a straight 
line passing through the points (pi,po) and {pf,pQ + e), 
where we see that for a small e [e — {pf — Pi) tan(a)] a 
good approximation of the line Pr = Po — constant is 
found. We choose arbitrary values for the energy density 
interval {pi,Pf) = (1,2) and pressure po = 1, with these 
values the equation of state in the coexistence region is 
of the form, 

Pr = l + e{pr - 1) 1 < Pr < 2 . (21) 

When e ^ we obtain the case of constant pressure 
= 1. In figure 121 we solved numerically the TOV equa- 
tions (|20|l with the equation of state H21|l for different 
values of e. In the limit e ^ we can affirm that the 
transition occurs only in a concentric spherical surface 
of the star. This means that in the limit of e ^ for 
a given equation of state of the type H16|l with l|21|) . the 
TOV system brings us a solution of a star model that 
consists in a foliation of spherical, isobaric, and acronal 
surfaces. This is in accordance with the stationary hy- 
pothesis made for the model, because if the transition 
take place in a radial interval instead of a fixed radius, 
the different gravitational forces at different radius will 
make a non stationary system. 

IV. SOLUTION OF TOV EQUATIONS FOR A 
SYSTEM IN PHASE TRANSITION 

A. Thermodynamic properties of the star 

For didactic purposes we choose the isothermal with re- 
duced temperature Tr = 0.9 and parameters (x = 0.005, 
^ — 0.00005) to study the thermodynamic properties and 
the space-time behavior. It can be verified that all the 
isotherms that allow a phase transition have the same 
qualitative behavior. With the values of the coexistence 
pressure po — 0.583 and the energy density coexistence 
interval (pi, ps) = (0.350, 1.623) found from ifTH) . we can 
solve numerically the system of equations H20|) . The in- 
ner boundary condition can be set recalling that m(r) is 
proportional to the mass inside the sphere of radius r, 
m(r = 0) = 0. So, the inner boundary condition for the 

new variables are (r,po,Fo) = (0,2,0); which represent 



a liquid core of high energy density, the value po = 2 is set 
arbitrarily. In figure |21 we show the numerical solution 
of the TOV equations for the thermodynamic quantities 
in the star with phase transition. The phase transition 
shows up like a suddenly brake in the energy density at 
r « 0.6A. Looking at the energy density plot we can 
think the high dense region with r < 0.6A as a liquid 
core and the region with small density r > 0.6A as a gas 
mantle. The pressure has a discontinuity in its radial 
derivative at r « 0.6A and presents the expected behav- 
ior, i.e. high and low pressures appears in the high energy 
density (core) an low energy density (mantle) regions re- 
spectively. The geometric mass increases rapidly in the 
liquid region and slowly in the gas region, it also presents 
a discontinuity in its radial derivative at the same radius 
of the star. We note that the energy density and the 
pressure goes to zero for large values of r, so in principle 
we can match continuously to a Schwarzschild free-space 
solution, as was done in neutron-star models 16], or to 
an exterior dust fiuid (with zero pressure). 



B. Space-time geometry in a star in phase 
transition 

The metric components in the star are found through 
the functions {iy,X) from equations © and (|10|) . The 
integration constant from equation (|10|l can be obtained 
by imposing that the component 500 Sii a distant radius 
has to be equal to the Schwarzschild exterior solution 

5oo(r.) = l-^^, (22) 
rs 

where is the radius of the star. The component gn is 
calculated directly from equation l^. 

These components are depicted in figure 21 We note 
in this figure that the component goo of the metric does 
not present discontinuities in the phase transition region 
as found in the radial derivative of the gn component at 
r w 0.6A. This motivate us to think that, in general, the 
phase transition could be manifested only in the spatial 
part of the metric. 



C. Curvature associated with the phase transition 

To study the space-time curvature in a system with 
phase transition, we used the Kretschmann invariant 
{K = R°'P^^Ra(}^s). With the help of relations Q, we 
find that the Kretschamnn invariant can be written as 
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In figure 21 "we see that the Kretschmann invariant 
present a similar behavior than the energy density, i.e. 
highest curvature values in the core of the star and at long 
distances the curvature tends to zero. The Kretschmann 
invariant depends on the values of the pressure, energy 
density and geometric mass, so it is normal to expect a 
discontinuity in the coexistence region at the same sphere 
radius (r « 0.6A). For a better understanding of the cur- 
vature associated to the phase transition, we project the 
Riemann tensor along the orthonormal tetrad ^ , where 
a is the contravariant tensor index and (a) a label dis- 
tinguishing the particular vector. So, we have 

R{a){b){c){d) = -Ra/375CfQ)Cfb)?('c)Cfd), (24) 

in our case the orthonormal tetrad is 



state that the abrupt change in the Kretschmann scalar 
can be thought as a curvature transition in the star due 
to the loss of energy needed for the phase transition to oc- 
cur, see the energy density profile from figureOl We note 
that discontinuities in the scalar curvature are not rare 
phenomena in nature. One simple example of this dis- 
continuity is provided by a simple star model described 
in its exterior by the Schwarzschild solution and its in- 
terior by the solutions corresponding to a homogeneous 
sphere [121 . In this case there exist a discontinuity in the 
scalar curvature mainly because the energy 

density has a discontinuity, constant inside the star 
and zero outside. 



^fo) = (e-^/^ 0,0,0), 

^fi) = (0,6-^/2,0,0), 

Cf2) = (0,0,r-i,0), 

Cf3) = (0,0,0, [rsin(0)]-i). 



(25) 



and we obtain the following non null components of the 
Riemann tensor: 



^(o)(i)(o)(i) " A2 1 + y ~ ^ I ' 



^(0)(3)(0)(3) = ^(0)(2)(0)(2) 



^(1)(3)(1)(3) = ^(1)(2)(1)(2 




R, 



(2)(3)(2)(3) 




(26) 



These components are plotted in figure^ We see that all 
the components are functions of the class . We have 
discontinuous radial derivative at r w 0.6A. In particu- 
lar the components i?(i)(2)(i)(2) and i?(o)(i)(o)(i) present 
the same abrupt change that the Kretschmann scalar. 
The Kretschmann scalar can be written as a function of 
the non zero components of the Riemann-Christoffel cur- 
vature tensor So, we can assure that the abrupt 
change in the Kretschmann scalar is due to the com- 
ponents -R(i)(2)(i)(2) and i?(o)(i)(o)(i). These two compo- 
nents depend on the reduce energy density pr- So, we can 



CONCLUSIONS 



We studied a model of a stationary spherical symmet- 
ric system in phase transition. The system is a ball of 
fluid made with identical and auto- gravitating particles 
with the energy-momentum tensor of a perfect fluid that 
obeys a van der Walls like equation of state. This model 
appears to be the most simple and general model to de- 
scribe a stellar object that undergoes a phase transition. 
The phase transition in the thermodynamic quantities is 
manifested as an abrupt change in the energy density, 
and discontinuities in the radial derivatives of the pres- 
sure and geometric mass at a fix radius. This radius 
divides the star in two regions: one with a high energy 
density (a liquid core) and the other with low energy 
density (a gas mantle). In the space-time geometry of 
the star, we see an abrupt change in the Kretschmann 
invariant and some components of the Riemann tensor. 
These can be thought as a curvature transition process 
due to the loss of energy needed for the phase transition 
to occur. A discontinuity in the radial derivative is also 
present in the spatial component of the metric but not 
in its time component, this lead us to think that phase 
transitions could be manifested only in the spatial part 
of the metric. For the results obtained in this work, we 
interpret the stationary system studied as a foliation of 
concentric and isobaric spherical surfaces where the phase 
transition occurs in only one of these surfaces. We note 
that the procedure used in this article can be applied 
to other models, such as systems with different kind of 
particles and systems characterized with more than one 
phase transition. 
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FIG. 1: A typical van dcr Walls isothermal for a system in gas- liquid phase transition using reduce quantities. The coexistence 
pressure po and the coexistence region pi < pr < ps is calculated using the Maxwell construction. 





FIG. 2: Analysis of the energy density and pressure profiles. We see from the graphs that in the limit when e —^ {pr = 
constant), the energy density profile tends to a vertical line at a fixed radius while the radial domain of the pressure profile 
becomes a point. 




FIG. 3: Energy density, pressure and geometric mass profiles for the interior of a star in phase transition. The discontinuities 
in their radial derivatives at f « 0.6 represent the change of phases from a liquid core (high energy density) to a gas mantle 
(low energy density). 



0,1 0,2 0,3 0,4 OS 0.6 0,7 0.8 0,9 1 

Radius f 



0.1 0.2 0.3 0,4 0.5 0.6 0.7 08 0.9 1 

Radius r 




FIG. 4: Behavior of the goo and gn components of the metric, the Kretschmann curvature invariant and the components of 
the Riemann-Christoffel tensor in a star in phase transition. We note that the phase transition is manifested only in the spatial 
component of the metric. The Kretschmann invariant, and curvature components -R(i)(2)(i)(2) and i?(o){i){o){i) experience a 
sudden change, this can be interpreted as a curvature transition inside the star due to the loss of energy needed for the phase 
transition to occur. 



